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Abstract 

We calculate particle spectra and continuum photon emission from the Cassiopeia A supernova rem- 
0> , nant (SNR). The particle spectra, ion and electron, result from diffusive shock acceleration at the 

^^ ' forward SNR shock and are determined with a nonlinear Monte Carlo calculation. The calculation 

self-consistently determines the shock structure under the influence of ion pressure, and includes a 
simple parameterized treatment of electron injection and acceleration. Our results are compared to 
photon observations, concentrating on the connection between the Radio and GeV-TeV gamma-ray 
range, and to cosmic ray ion observations. We include new upper limits from the Cherenkov Array 
C^^ . at Themis (CAT) imaging Cherenkov telescope and the Whipple 10m 7-ray telescope at > 400 GeV. 

These new limits support the suggestion (e.g., Cowsik & Sarkar 1980; Allen et al. 1997) that ener- 
^ ' getic electrons are emitting synchrotron radiation in an extremely high magnetic field (~ 1000/uG), 

far greater than values routinely assigned to the ISM, and help to constrain our model. The large 
magnetic field allows acceleration of cosmic ray ions to well above 10^^ eV per nucleon in the ~ 300 
yr lifetime of Cas A. 
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Q\ '. 1 Introduction 

We describe particle acceleration in an expanding, spherical SNR blast wave with a plane-wave, 
Qh| steady-state shock model. The justification for using a steady-state calculation to model time- 

Q ■ dependent SNRs is given in Ellison & Berezhko (1999) and the full details and assumptions of our 

model are given in Baring et al. (1999). Briefiy, the global SNR shock parameters (e.g., shock speed 
and radius as a function of remnant age) used as input for the model, are estimated with a simple 
Sedov solution of the evolving blast wave in an uniform medium. Given the global shock dynamics, 
we are able to calculate the shock structure and first-order particle acceleration assuming that elec- 
trons and ions are accelerated directly by the forward shock, leaving considerations of the reverse 
C^ ' shock, where substantial X-ray emission may originate, for later work. Important limitations of our 

current model is that we do not consider oblique magnetic field structures or include second-order 
Fermi acceleration, which may be important in the low Alfven Mach number shocks {Ma_i ~ 4 for 
the parameters we show here) implied by the large i?-fields (see Bykov & Uvarov 1999 for a model 
of electron injection which does include 2nd-order acceleration). 

For the particular case of Cas A, the forward shock may be interacting with pre-supernova wind 
material swept into a relatively dense shell (e.g., Borkowski et al. 1996) which may account for 
the high magnetic field if the stellar magnetic field was compressed along with the wind material. 
Alternatively, the large field may be the result of amplification by turbulent eddies (Keohane 1998 
and references therein). If magnetic fields B r^ 1000/iG dominate the acceleration region, ions can be 
accelerated to well above 10^^ eV/nuc. While the overall fiuxes of these particles may not be sufficient 
to provide the bulk of the cosmic rays at 1-10 GeV because of the small numbers of particles swept up, 
particle spectra from young SNRs might be harder than from older, slower SNRs and may dominate 
above ~ 10^'* eV through the 'knee' near 10^^ eV. 
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Figure 1: Cassiopeia A spectrum. Observations are from Baars et al. (1977: radio); Lagage et 
al. (1999, in preparation: ISOCAM); The et al. (1996: RXTE and OSSE); Esposito et al. (1996: 
EGRET); Lessard et al. (1999: Whipple); and Goret et al. (1999: CAT). The model photons come 
from a single set of proton, helium, and electron spectra calculated with unshocked parameters shown 
in the figure. A single normalization has been applied to all components to match the radio flux. The 
subscript '1 ' indicates unshocked values and Msi and Mai o-tc the sonic and Alfven Mach numbers. 
Note that the bremsstrahlung emission cuts off at a much lower energy than the pion-decay due to 
the synchrotron losses the electrons experience. In these preliminary results, the inverse Compton 
does not included synchrotron- self- Compton although this may dominate the inverse Compton from 
the primordial background radiation. 

Not only does the high S-field make it possible to produce cosmic rays to > 10^^ eV, but a 
homogeneous model with a single set of parameters can be found which affords a reasonable fit 
to the intensities of diffuse photon emission from radio to 7-rays. Limits imposed by the radio 
and 7-ray observations allow us to place constrains on the electron-to-proton ratios produced by 
shock acceleration. The detailed shapes of individual components (e.g., radio, X-ray), however, are 
difficult to model with a single set of parameters. Despite this limitation, important model and/or 
environmental constraints can be inferred if it is assumed that the relativistic electrons responsible 
for the radio synchrotron emission also produce the diffuse infrared and X-ray continuum via the 
synchrotron mechanism, (and in the same emitting volume). Likewise, if the GeV and TeV 7-ray 
emission is dominated by inverse Compton and bremsstrahlung emission (rather than pion-decay 
from energetic ions), it can be assumed that these same relativistic electrons are responsible for the 
GeV and TeV gamma rays. Assuming a homogeneous environment precludes, of course, the modeling 
of emission from the lumpy morphology, knots, ring structure, etc., that characterizes high spatial 
resolution observations of Cas A. 

2 Results 

In Fig. 1 we compare our results to observations of Cas A from radio to TeV 7-rays. We have a 
single normalization for all of the model components and this has been chosen to match the radio flux. 
Since the same electron distribution that produces the radio emission can also produce bremsstrahlung 
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Figure 2: Cosmic ray spectra measured at the Earth (data from, Shibata 1996, 1998) compared to 
model spectra with sam,e param,eters as shown in Fig. 1. The model spectra have been multiplied by 
Rigidity^^-^^ to account for escape during propagation and have been normalized separately to match 
the observations. The light dotted histogram is the model Fe spectrum without nuclear destruction 
during propagation while the dot-dashed histogram includes nuclear destruction. The turnover in the 
data below ~ 1 Ge V/A is from solar modulation, which is not included in the model. 

and inverse Compton emission at GeV-TeV energies, the combination of the observed radio intensity 
and the low 7-ray upper hmits forces the conclusion that the magnetic field is exceedingly high. 
Otherwise, the bremsstrahlung and inverse Compton emission at TeV energies would have been 
observed. This might not be the case if the radio emitting electrons occupy a greater emission 
volume than GeV-TeV emitting electrons, or that the lower energy (i.e., radio-emitting) electrons 
can preferentially sample regions of clumped magnetic field and/or density. Even if either of these 
cases arises, the inverse Compton component, which depends at most only weakly on the background 
particle density, will set a lower limit on B that remains well above standard ISM values. Note that 
the cosmic microwave background radiation forms the seed for inverse Compton scattering in the 
present exposition; work is in progress to include synchrotron self-Compton contributions, which will 
tighten the constraints discussed here, implying even higher S-fields. 

The relative importance of the pion-decay emission depends sensitively on the model parameters 
and we have not yet done a careful survey of the parameter space. For our preliminary results shown 
in Fig. 1, we have chosen our electron injection parameters to give an electron to proton ratio at 
relativistic energies, (e/p)ioGeV — 0.08, somewhat above observed cosmic ray values. 

We again emphasize that even though our model is for plane-parallel, steady-state shocks, de- 



tailed comparisons (Ellison &: Berezhko 1999) with the spherically symmetric, time-dependent model 
of Berezhko (1996; see also Berezhko, et al. 1996) show that the steady-state and plane shock approx- 
imations do not seriously influence the results as long as the diffusion length of the highest energy 
particles is a small fraction of the shock radius, i?sk, as should be the case in Cas A and other young 
SNRs. If Cas A is currently interacting with a relatively dense shell of material formed from the 
pre-SN wind (Borkowski et al. 1996), the Sedov solution can be replaced by estimates of the ambient 
density, ni, magnetic field, Bi, i?sk, and Vg]^, which translate to maximum particle energies. 

3 Conclusions 

The observed radio intensity of Cas A, combined with the EGRET and TeV upper limits, imply 
magnetic fields ^ 1000//G. Fields this large make it possible to accelerate cosmic ray ions to above 
10^^ eV/nuc in the ~ 300 yr lifetime of the remnant, since the time to shock accelerate ions of charge 
Q to -Emax is (e.g., Baring et al. 1999) 

iacc ^ 190 i^-j (2000km/sj [W^J [w^J ^' ' 

Here, r] is the number of gyroradii in a scattering mean free path and is approximately one in the 
Bohm limit, which we assume in this model. In Fig. 2 we compare the model spectra, all multiplied 
by i?""'^^ [R = pc/{Qe)] to model rigidity-dependent escape during propagation, with cosmic ray 
observations. No attempt has been made to model the abundances of these components, each being 
separately normalized to match the observations. It's clear that the spectra extend through the knee 
region. Young SNR shocks sweep up far less ISM material than older, larger shocks, but the high 
energy cosmic rays produced may have flatter spectra than the bulk of the cosmic rays accelerated by 
older SNRs which have weaker shocks and, if so, could dominate the cosmic ray flux near the knee. If 
this is the case, spectral (and perhaps compositional) features should exist in the cosmic ray spectra 
as these components become dominant. 

If high i?-fields are common in young SNRs, it has important implications for 7-ray emissivity as 
well as cosmic ray production. High fields imply that radio intensities will be high, concomitant with 
relatively low relativistic electron fluxes (and presumably low relativistic ion fluxes as well), lowering 
the 7-ray emissivity. Therefore, radio loud SNRs may not be the best candidates for 7-ray studies, 
and some other indicator may be required to guide observational programs. 

Allen, G.E. et al. 1997, Ap. J. (Letts), 487, L97 

Baars, J.W.M., Genzel, R., Pauliny-Toth, I.I.K., & Witzel, A. 1977, A&A, 61, 99 

Baring, M. C, Ellison, D. C, Reynolds, S. R, Grenier, I. A., & Goret P. 1999, Ap.J., 513, 311 

Berezhko, E.G. 1996, Astroparticle Phys., 5, 367 

Berezhko, E.G., Yelshin, V.K. & Ksenofontov, L.T. 1996, ZhETF, 109, 3. 

Borkowski, K.J., Szymkowiak, A.E., Blondin, J.M., &: Sarazin, C.L. 1996, Ap.J., 466, 866 

Bykov, A.M., & Uvarov, Yu.A. 1999, JETP, 88, 465 

Cowsik, R., & Sarkar, S. 1980, M.N.R.A.S., 191, 855 

Ellison, D.C., & Berezhko, E.G. 1999, 26th ICRC (Salt Lake City), OG 3.3.27. 

Esposito, J. A., Hunter, S.D., Kanbach, G., &: Sreekumar, P. 1996, Ap.J., 461, 820 

Goret, P. et al. 1999, 26th ICRC (Salt Lake City), OG 2.2.18. 

Keohane, J.W. 1998, Ph.D. Thesis, University of Minnesota. 

Lessard, R. (Whipple Collaboration), 1999, Proc. 19th Texas Symposium, Paris 1998, in press. 

Shibata, T. 1996, Nuovo Cimento C, 19, 713. 

Shibata, T. 1998, private communication. 

The, L.-S., et al. 1996, A. A. SuppL, 120, 357 



